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Abstract. We have used global magnetohydrodynamic simulations of the interaction 
between the solar wind and magnetosphere together with single particle trajectory 
calculations to investigate the sources of plasma entering the magnetosphere. In all of our 
calculations solar wind plasma primarily enters the magnetosphere when the field line on 
which it is converting reconnects. When the interplanetary magnetic field has a 
northward component the reconnection is in the polar cusp region. In the simulations 
plasma in the low latitude boundary layer (LLBL) can be on either open or closed field 
lines. Open field lines occur when the high latitude reconnection occurs in only one cusp. 
In the MHD calculations the ionosphere does not contribute significantly to the LLBL for 
northward IMF. The particle trajectory calculations show that ions preferentially enter in 
the cusp region where they can be accelerated by non-adiabatic motion across the high 
latitude electric field. For southward IMF in the MHD simulations the plasma in the 
middle and inner magnetosphere comes from the inner (ionospheric) boundary of the 
simulation. Solar wind plasma on open field lines is confined to high latitudes and exits 
the tailward boundary of the simulation without reaching the plasma sheet. The LLBL is 
populated by both ionospheric and solar wind plasma. When the particle trajectories are 
included solar wind ions can enter the middle magnetosphere. We have used both the 
MHD simulations and the particle calculations to estimate source rates for the 
magnetosphere which are consistent with those inferred from observations. 

1. Introduction 

There are two main sources of magnetospheric plasma: the solar wind and the 
ionosphere. Classical theory of the solar wind interaction with the magnetosphere 
predicts that the magnetopause should be an impenetrable boundary separating the 
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shocked solar wind plasma of the magnetosheath from the hot tenuous plasma of the 
magnetosphere. However, three decades of observations have demonstrated clearly that 
magnetosheath plasma exists inside all regions of the magnetopause (see Sibeck et al. 
[1999] for a recent review). These magnetospheric boundary layers include the low- 
latitude boundary layer (LLBL), the polar cusp entry layer and the high latitude plasma 
mantle (PM). 

Observations from low-altitude spacecraft indicate that shocked solar wind 
plasma from the magnetosheath enters the magnetosphere over a wide region under 
almost all solar wind conditions [e.g. Newell and Meng, 1992]. Observations in the 
magnetotail lobes indicate that plasma enters the magnetosphere along its entire length. 
Recently Hultqvist et al. [1999] have combined the results from a large number of studies 
to identify and quantify the major source regions for magnetospheric plasma. They 
estimate that ~10 26 s' 1 solar wind particles enter through the dayside magnetopause while 
10 -10 s enters along the tail magnetopause. However most of the particles that enter 
through the tail magnetopause also exit via the tail and never reach the closed field line 
region of the magnetosphere where the required source strength is only a few times 
10 /s. Similarly, Moore et al. [1999] organized the published observations from polar 
orbiting spacecraft [Yau et al., 1985; Abe et al., 1996] to summarize the ionospheric 
outflow rates. Their results show that the ionosphere can provide ~10 26 /s of 0 + and H^ to 
the magnetosphere. 

A number of processes have been proposed to account for the transfer of solar 
wind plasma into the magnetosphere. These include magnetic reconnection, finite Larmor 
radius effects, diffusion, the Kelvin-Helmholtz instability, impulsive penetration and 
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direct cusp entry. In this paper we use computer simulations to study the sources of 
plasma for the magnetosphere with emphasis on the sources of the LLBL. First we 
investigate the sources of plasma by using a global magnetohydrodynamic (MHD) 
simulation of the interaction of the solar wind with the magnetosphere. We use the MHD 
model to investigate both the ionospheric and solar wind sources. Then we add single 
particle trajectory calculations of solar wind ions to the global MHD calculations in order 
to estimate the importance of phenomena such as finite Larmor radius effects not 
included in MHD. For each approach we estimate the source rate. In section 2 we briefly 
describe the models and techniques used in this study. The simulation results for 
northward IMF are presented in section 3 while those for southward IMF are in section 4. 
Finally in section 5 we discuss these results in the context of magnetospheric 
observations. 

2. Approach 

Since the details of the global magnetohydrodynamic simulation code used for 
this study have been described previously [Ogino et al., 1992; 1994] we will describe it 
only briefly. We solve the resistive MHD equations and Maxwell’s equations as an initial 
value problem on a rectangular grid. For the cases with purely northward and southward 
IMF we used a (402x102x102) points grid with a uniform mesh size of 0.4R E . For these 
calculations the simulation box extends from 30 R E >x> -130R £ , 0 < y < -A0R E , and 
0 <z < 40 R e . We also have included results from simulations with IMF B Y ^ 0 . For 
these cases we used a simulation box with 322x82x162 points and a mesh size of 0.5R E . 
In the Z direction the calculation domain is - 40 R B <z< 40 R E . We solve the 
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differential equations by using a modified version of the Leapfrog scheme that is a 
combination of the Leapfrog scheme and the two-step Lax-Wendroff scheme. The 
simulation parameters are fixed to solar wind values at the upstream edge of the 
simulation box with free boundary conditions at the sides and back. For all the 
simulations we used symmetry boundary conditions at y = 0 while for the quarter 
magnetosphere calculation we also used symmetry boundary conditions at z-0. The 
ionospheric boundary is at 3.5R E . The solar wind velocity was 300km/s, the density was 
5cm' 3 and the temperature was 2x1 0 5o K. 

The simulation for northward IMF was initialized by using an unmagnetized solar 
wind flow that lasted for up to three hours. Then a northward IMF was introduced into 
the flow. After four hours with northward IMF during which time a quasi-steady 
magnetosphere developed, the IMF was turned southward and the simulation was run for 
another four hours. For the simulations with By^O the initialization period was followed 
by a five-hour interval with northward IMF. Then the IMF was rotated in 15° increments 
in the YZ plane with a dwell time at each direction of 15 minutes [ Walker et al., 1999]. 
While the dwell time was sufficient to enable us to determine the entry mechanism, the 
overall magnetospheric configuration contains the effects of changing IMF orientations 
[Walker et al., 1999], For selected IMF directions the dwell time was increased to 2 hours 
to enable us to investigate the configuration of a quasi-steady magnetosphere. For all of 
the simulations B = 5nT . 

We also examined the entry of solar wind ions into the magnetosphere by 
calculating the ion trajectories using a fixed magnetic and electric field model [Ashour- 
Abdalla et al., 1993], We use magnetic and electric fields from the MHD simulations. 
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The electric field is given by E = — v xB + r/J where v is the velocity, B the magnetic 

field, and J is the cinrent density. An explicit resistivity (r| ) is included in the MHD 
simulations [ Ogino et al., 1994], In each case a distribution of ions with the temperature 
used in the MHD simulations was launched in the solar wind upstream of the bow shock. 
The ion trajectories were calculated by solving the equation of motion by using a fourth 
order Runge-Kutta method. 

3. Northward IMF Results 

We started our studies by modeling plasma entry when the interplanetary 
magnetic field was northward. In Figure 1 we have plotted flow streamlines for a case 
with purely northward IMF. We show one quarter of the magnetosphere with pressure 
contours on the noon-midnight meridian and equatorial planes. The bottom panel 
contains a blow up of the region near the Earth. This simulation was run for 3 hours with 
no IMF and then for 4 hours with a northward IMF. This snapshot was taken 2 hours after 
the northward IMF entered the simulation box and after a quasi-steady magnetospheric 
configuration had developed [Ogino et al., 1994; Bargatze et al, 1999], Most of the flow 
impinging on the magnetosphere is diverted around the obstacle (top). However, some 
plasma can enter the magnetosphere. In the bottom panel we have plotted additional 
streamlines from the region around the subsolar point. These streamlines enter the 
magnetosphere and move tailward on magnetospheric field lines. In Figure 2 we have 
reproduced two of these flow streamlines and calculated magnetic field lines along each 
streamline. The red field lines correspond to a streamline along the noon-midnight 
meridian slightly northward of the equator. This streamline moves from solar wind to 
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magnetospheric field lines at location 3 when the field line on which it is convert ing 
reconnects in the polar cusp region forming a new closed magnetospheric field line. The 
streamline then moves around the flank magnetosphere into the near-Earth tail region. 
Note that the streamline remains on field lines that map to the LLBL while on the 
dayside. The second streamline starts along the noon-midnight meridian at the equator. 
The corresponding field lines are blue. Again the streamline enters the magnetosphere 
when the field line on which it is converting reconnects with high latitude tail field lines. 
The streamline then moves along the dayside LLBL and into the tail. 

The magnetospheric configuration changes significantly when the IMF has a By 
component. In the MHD model for purely northward IMF, the IMF field lines reconnect 
simultaneously in the northern and southern hemisphere forming closed field lines. When 
there is a By component the reconnection can occur in only one hemisphere. In Figure 3 
we have plotted pressure contours in the equatorial plane along with selected magnetic 
field lines from a snapshot from a simulation when the IMF had northward and 
dawnward components. At this time the IMF had been northward for 45 minutes and 
duskward for 90 minutes. The pressure contours have been made transparent so that field 
lines (in purple) south of the equator can be seen. Note the open field lines shown here 
along the dusk flanks. A flow streamline starting in the southern hemisphere just 
duskward of noon has been drawn in orange. As in the purely northward IMF case, this 
streamline enters the magnetosphere at approximately position 2 when the field line on 
which it is convecting reconnects at high latitude. The streamline enters the 
magnetosphere on open field lines and remains on open field lines until it exits the 
simulation box. 
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We also launched a distribution of solar wind test particles in the magnetic and 
electric fields from this MHD simulation. The distribution function of the test particles 
was selected to have the same plasma moments as the MHD simulation. The test particles 
were all protons and were launched in the plane at jc = 20 R E . The white trace in Figure 3 
gives the trajectory of one of these ions. This proton has energy of about 190eV in the 
solar wind frame. The ion entered the magnetosphere at approximately position 2. Again 
the entry mechanism is high latitude reconnection. Using similar calculations Richard et 
al., [1994] determined that high latitude reconnection is the main way in which ions enter 
the magnetosphere when the IMF is northward. After entering the magnetosphere the ion 
bounced several times at high southerly latitudes before moving onto open field lines and 
then into the polar cusp region. The ion emerged from the cusp onto closed dayside field 
lines. In Figure 4 this trajectory has been plotted from two perspectives. The view in the 
top panel is from the afternoon while the view in the bottom is from the morning. Field 
lines along the trajectoiy have been color coded to indicate the magnetic field magnitude. 
Note that soon after the ion crossed the magnetopause it bounced about the off equatorial 

minimum B that is characteristic of the magnetosphere just earthward of the dayside 
magnetopause. This type of oscillation about the high-latitude minimum B has been 

reported by Delcourt et al. [1992] based on particle calculations in empirical E and 

B field models (see their Figure 4). However, by the time the ion in Figure 4 reached the 
morning side (bottom) it was on a trapped orbit bouncing about the equator. In between 
the ion moved straight across the dayside cusp near the region where the high latitude 
reconnection is occurring. 
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The effects of the motion across the cusp can be seen in Figure 5. Here the 
trajectory has been projected onto three orthogonal planes. Each trajectory is color coded 
to show the energy of the particle. While the ion is on the dusk side it has about the same 
energy as it had in the magnetosheath. However, when it moved across the polar cusp it 
quickly gained ~20kV. Richard et al., [1994] argued that ions move non-adiabatically in 
the polar cusp region when high latitude reconnection causes the radius of curvature of 
field lines there to become comparable with the ion gyro-radius. This non-adiabatic 
motion carries the ions across the high latitude electric field where they are accelerated. 

The ion in Figure 5 is unusual in two ways. First it gained more energy in the cusp 
region than typical ions. Second it entered the cusp region after moving in the dusk side 
magnetosphere. However, our trajectory calculations indicate that non-adiabatic motion 
across the polar cusp region is the main way in which the most energetic ions are 
accelerated. Most of the accelerated ions enter the cusp directly from the magnetosheath 
when the field line on which they are moving reconnects at hi g h latitudes. 

In Figure 6 we have plotted the trajectory of an ion that did not enter the high 
latitude reconnection region. This particle entered the magnetosphere very close to noon 
when the field line on which it was moving reconnected with the southern lobe field (3). 

It quickly moved to closed field lines and bounced and drifted near the dayside 
magnetopause to the night side (4-6). In Figure 7 we have plotted the energy of this 
particle along its trajectory. This ion gained less than one keV as it moved across the 
dayside magnetosphere. 
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Table 1 


Ion Entry Rates 2?z>0 



MHD 

Ion Trajectory 

By~ 0 

4x10 26 /s 

lxl0 27 /s 

By^O 

5x10 26 /s , 5x10 27 /s 

2x10 27 /s, 5x10 27 /s 


In Table 1 we have estimated the ion entry rates for the northward IMF 
simulations. We estimated the entry rate in the MHD simulation by determining the 
streamlines that cross the magnetopause onto magnetospheric field lines and 

calculating j/? v-ds where p is the mass density, v is the fluid velocity and ds is an area 

element on the magnetopause. The integral was taken over the magnetopause area 
threaded by penetrating streamlines. Since the variable in the MHD simulations is the 
mass density, we assumed the solar wind to consist solely of protons in order to obtain 
the entry rate. When the IMF is purely northward 4X10 26 /s enter on closed field lines. 
When the IMF has a By component we have listed two numbers. The first is the entry rate 
on closed field lines. It remains about the same as in the case with B YIMF = 0 . The second 

number is the entry rate on both open and closed field lines. It is an order of magnitude 
larger because of ions that are inside the magnetopause but on open field lines. For 
th eB Y * 0 case we used a simulation with B r = -3.5 nT and B Y = 3.5nrthat was run 
until a quasi-steady magnetosphere was formed (2 hours). 

We also estimated the entry rate by using the particle trajectory calculations. The 

entry rate is given by r = x Vsiv x Hs,v ^y where C is the number of ions entering the 
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magnetosphere, ATis the total number launched, A is the area in the solar wind over which 
the ions were launched, vsw is the solar wind velocity and nsw is the solar wind number 
density. In general the source rate estimates are larger from the trajectory calculations. In 
Figure 8 we have plotted the location on the magnetopause where the ions entered the 
magnetosphere. The entry points are bunched near the polar cusp. In MHD the 
comparable distribution has uniform entry. We get a higher estimate of the source rate 
from the particle calculations because of the enhanced entry in the cusp that is caused by 
the non-adiabatic motion [. Richard et al, 1994], 

4. Southward IMF results 

We have plotted flow streamlines and pressure contours for a case with purely 
southward IMF in Figure 9. This snapshot was taken 60 minutes after the IMF was turned 
southward from northward. By this time subsolar point reconnection had been occurring 
for approximately 45 minutes [Walker et al., 1993; Bargatze et al., 1999], The 
magnetosphere had not reached a steady state. A few minutes (t=70 minutes) after this 
snapshot reconnection began in the near-Earth (x « -15i? £ ) magnetotail. About 10 
minutes later reconnection began on lobe field lines. Since the magnetospheric 
configuration was not steady we calculated flow streamlines through a series of time 
steps starting before and extending after this time. Since the results were very similar to 
those in Figure 9 we have not included them in this paper. 

In the simulation flow streamlines from the solar wind do not reach the inner parts 
of the magnetosphere. The streamlines that reach the middle and inner parts of the 
magnetosphere come from the inner boundary of the simulation. Streamlines in the LLBL 
come from both the solar wind and the inner boundary. Those nearest the Earth come 
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from the inner boundary while streamlines nearer the magnetopause come from both 
sources. Note that the earthward flowing streamlines from the inner boundary turn around 
in the dayside magnetosphere and join with streamlines from the solar wind in flowing 
tailward. This can be seen more clearly in Figure 10. Here we have plotted only three 
streamlines along with magnetic field lines along the streamlines. The two sources of the 
LLBL have blue and green streamlines. Plasma on the blue streamlines enters the 
magnetosphere when the field line on which it is convecting reconnects at local times 
away from noon (2). This plasma moves on open field lines through the LLBL into the 
tail (3-5). In this example plasma from the inner boundary has taken a less direct path to 
the LLBL. The streamline emerges from the inner boundary on the dayside (green 1) and 
moves toward the dayside magnetopause. There it is picked up by reconnection (2) and 
moves on open field lines into the tail (3-5). In the tail it is picked up by earthward 
convection (6) and moves back toward the dayside magnetosphere (6-8). At the dayside 
magnetopause (between 8 and 9) the streamline again encounters reconnection. Finally it 
moves tailward on open field lines (10). In this example we followed a flow streamline 
from the inner boundary on the dayside. However plasma also emerges from the inner 
boundary on the nightside (Figure 9). These streamlines come from auroral latitudes. 
When the nightside flow streamlines reach the equator they too flow earthward into the 
dayside magnetosphere. Most of the outflow at the inner boundary is field aligned. 

The red field lines in Figure 10 have been calculated along a streamline that 
originates in the solar wind along the Earth-Sun line. It too enters the magnetosphere 
because of dayside reconnection (2) but then moves on open field lines the length of the 
simulation box (3-5). 
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We also launched a distribution of ions from the solar wind into the electric and 
magnetic fields from this simulation. We have plotted an example showing the trajectory 
of one of those particles in Figure 1 1 . This ion was launched along the Earth-Sun line and 
had an initial energy of 17eV in the solar wind frame. This ion had a small southward V| . 

Even in the solar wind the ion diverged slightly from the corresponding flow streamline. 
This became much more pronounced after it crossed the bow shock. The ion entered the 
magnetosphere when the field line on which it was moving reconnected (2). The ion 
moved along open field lines into the magnetotail (3-5) and reached the equator at 
x » —50 R e (6). Then it drifted earthward (7) and around the Earth before exiting the 
magnetosphere. By location (8) the field lines on which the ion was moving were in the 
LLBL very near the magnetopause. Finally the ion left the magnetosphere when the field 
line on which it was moving reconnected. 


Table 2 

Ion Entry Rates Bz<0 



MHD 

Ion Trajectory 

Solar Wind- 

X 

o 

& 

to 

2xl0 28 Is 

Open Field Lines 



Solar Wind - 


2xl0 27 Is 

Closed Field Lines 



Solar Wind - Direct Entiy 


clxlO 26 Is 

Ionosphere 

3X10 26 Is 
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We also estimated the ion source rate for this simulation. The results are 
summarized in Table 2. We calculated the solar wind source rate using the same method 
we used for the northward IMF case. In the MHD calculation about lxlO 28 /^ solar wind 
ions were entering the magnetosphere on open field lines at 60 minutes. As noted above 
this plasma exited the magnetotail boundary of the simulation without entering the closed 
field line region. That region of the magnetosphere was populated from the inner 
boundary of the simulation. We calculated the ionospheric source rate by integrating 

\p v-ds over a closed surface at the inner boundary. This gave 3 x 10 26 /s if we assume 

protons and 2 x 10 25 Is if 0 + . From the trajectory calculation we estimate that 2 x 10 28 /s 
solar wind ions enter on open field lines but that only about one in ten reach closed field 
lines. Finally a very few ions enter the magnetosphere by drifting through the 
magnetopause. Since the number of these particles is very low we could only obtain an 
upper limit of lxlO 26 Is for this source. 

5. Discussion 

In all of our calculations solar wind ions primarily enter the magnetosphere when 
the field line on which they are convecting reconnects at the dayside magnetopause. For 
northward IMF the reconnection is in the polar cusp and these particles form the low 
latitude boundary layer. For B z >0 and B Y *0 these ions can be on either open or closed 
field lines. Open field lines occur when the high latitude reconnection occurs in only one 
cusp. Since cases with simultaneous or nearly simultaneous reconnection in both cusps 
are thought to be rare [Fuselier et al., 2001], this is probably the way LLBL ions usually 
enter the magnetosphere. These ions quickly move from open to closed field lines. In the 
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MHD simulation the ionosphere does not contribute significantly to the LLBL for 
northward IMF. 

The situation is much more complex for southward IMF. In the MHD calculation 
solar wind plasma enters the magnetosphere on reconnecting field lines but never makes 
it to the closed field line region. Instead the open field lines convect to the end of the 
simulation box. Along the flanks of the magnetosphere open field lines carrying solar 
wind plasma are found in the same region as open field lines carrying ionospheric 
plasma. This region forms the outer part of the LLBL. Nearer the Earth the plasma in the 
MHD simulation is dominated by ionospheric plasma. When closed field lines carrying 
this plasma approach the magnetopause they can reconnect adding the ionospheric 
plasma to the LLBL mix. In the MHD simulation the inner regions of the LLBL are on 
closed field lines containing this ionospheric plasma. However, when we include the 
particle trajectory results we find that solar wind plasmas can reach the closed field 
region. Thus for southward IMF we would expect the LLBL to contain plasmas of both 
ionospheric and solar wind origin. 

Very recently Siscoe et al., [2001] have examined flow streamlines from an MHD 
simulation with Bz<0. They too found that solar wind plasma on reconnected field lines 
passes through the back of the simulation box (compare their Figure 5 with Figure 9 
(top)). In their simulation the plasma from the cusp region fills the plasma mantle and 
outer plasma sheet. They call the boundary between these two types of behavior the 
“fluopause” The “cusp” plasma in their simulation comes from the ionospheric boundary 
much like it does in our simulation (compare their Figure 7 with Figure 9 (bottom)). The 
“cusp” particles in both simulations populate the more distant parts of the tail. In our 
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calculation we continued to follow the streamlines throughout the magnetosphere and 
found that those in the middle tail return toward the Earth and end up populating the 
LLBL as well. We also found that the inner (ionospheric) boundary provides the plasma 
for the near magnetotail. This plasma comes from the nightside ionosphere boundary 
(Figure 9 top). Siscoe et al., [2001] did not present streamlines in this region. Overall the 
agreement between the two simulations is excellent. It is reassuring when two codes 
using very different approaches get the same results. 

Winglee [1998; 2000] has developed a multifluid MHD simulation that includes 
ionospheric outflow as a separate fluid. He too reports that the ionospheric plasma is 
confined to the inner magnetosphere for northward IMF while for southward IMF 
ionospheric plasma overcomes gravity and fills all but the more distant parts of the tail. 

Peroomian [2001] has mapped the locations where solar wind ions reach the 
magnetopause by calculating ion trajectories in a global MHD model. He used a 
somewhat stronger solar wind and found a few times 10 28 /s reached the magnetopause 
and that more ions reached the magnetopause for purely southward IMF than when there 
was an IMF By. He did not calculate the number entering the magnetosphere so we 
cannot compare the population of various regions of the magnetosphere. 

The source rates from the simulations in Tables 1 and 2 are in reasonable 
agreement with inferences from observations [Hultqvist et al., 1999]. When we use the 
trajectory calculations to estimate the solar wind source rates we get consistently higher 
source rates than with the MHD simulations. The largest difference occurs when the IMF 
is northward. This is primarily the result of entry into the polar cusp. Ions can undergo 
non-adiabatic acceleration in the cusp regions. In our simulations ions gained up to 
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20keV but acceleration to a few keV is more common. Finally we should remind the 
reader that the particle trajectory calculations are not self-consistent. The particles’ 
motion cannot influence the electric and magnetic fields through which they move. In 
reality the fields will change and so will the trajectories. Thus the trajectory plots shown 
in this paper should be thought of as a simple approximation to the actual motion. 

Moore et al., [1999] have summarized the ionospheric outflow to the 
magnetosphere based on observations from polar orbiting satellites (see their Table 2.1). 
They find outflows from lxl0 25 /s - lxl0 26 /s for FT and2xl0 25 /s - 2xl0 26 /s for 0 + with 
larger values for active times at solar maximum. The values in Table 2 fell at the low end 
of this range. The low energy (<100eV) observations used by Moore et al., [1999] come 
from low altitude (6000- 10,000km) spacecraft while the high energy observations (lOeV- 
17keV) are from high altitudes (16,000-24,000km). Since the low energy ions can be 
energized as they move upwards the observations may be over estimates. In our 
calculations the ionospheric source is passive. While we do not specify any ionospheric 
outflow the ionospheric boundary condition provides a cache of plasma that can be 
moved outward. The ionospheric outflow results solely from the solar wind 
magnetosphere interaction in our calculation. 

We have not carried out calculations in which we launch ionospheric ions in the 
global MHD fields. Other [Delcourt et al., 1992; 1994; Peroomian and Ashour-Abdalla, 
1996] have followed the trajectories of ionospheric ions in empirically based electric and 
magnetic field models. Delcourt et al., [1992; 1994] have followed protons from the 
tailward edge of the polar cusp into the magnetotail. The trajectories of their lowest 
energy particles (<100eV) are very much like the tailward part of the ionospheric 
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streamlines in Figure 9. Peroomian and Ashour-Abdalla [1996] launched IT, He + , and 0 + 
ions from the auroral zone. These ions populated the near-Earth tail region and the 
dayside magnetosphere. In the MHD model the auroral zone streamlines in Figure 9 (top) 
cover approximately the same region in the tail but don’t extend to the magnetopause on 
the dayside. 

Whether the IMF is northward or southward the simulation results indicate that 
reconnection is responsible for most plasma entry into the magnetosphere. While few 
would argue that reconnection is the dominate mechanism driving magnetospheric 
dynamics since its discovery [Hones et al., 1972; Eastman et al., 1976] it has been 
tempting to associate the LLBL with a viscous component of the solar wind 
magnetosphere interaction as well. Suggestions for the viscous stress include unstable 
Kelvin-Helmholtz (K-H) waves [Lee et al., 1981; Pu and Kivelson, 1983] that in 
principle can be studied with global MHD simulations. Two dimensional MHD models 
have been very successful in delimiting the flow and plasma conditions under which the 
K-H instability can grow on the dayside magnetopause [Muira, 1995] and along the 
flanks of the magnetotail [Otto and Fairfield, 2000], In the Otto and Fairfield [2000] 
study in the nonlinear stage of the K-H vortices, the magnetic field becomes highly 
twisted and local regions of oppositely directed field develop at which local reconnection 
allows plasma mixing across the magnetopause. To our knowledge only one 3D global 
MHD model has reported K-H waves. Walker et al., [1998] modeled the magnetosphere 
under conditions designed to maximize the chance of generating the K-H instability. 

They used a global MHD code to simulate the response of the magnetosphere to a density 
pulse (the density increased from 46 cm' 3 to 185 cm' 3 and back to 46 cm' 3 ) when the IMF 
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was strongly northward ( B z = 1 5nT ). Since the wavelengths of the K-H waves are on the 
order of the thickness of the LLBL (as small as 0 5R E ) a very fine grid was used 
( A* = 0.15^). The strong northward IMF was selected since Muira [1995] has shown 
that southward IMF stabilized the boundary and the density was increased to lower the 
Alfven velocity. After the pulse passed the magnetopause small K-H like flow vorticies 
were found along the flank magnetopause. We have examined the simulation in Figure 1 
for evidence of K-H waves but did not find any. Additional runs with higher resolution 
and varying solar wind parameters will be needed before we can use the simulation 
results to place limits on K-H waves as a source of LLBL plasma. 
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Figures 


Flow Streamlines and Pressure Contours 



Figure 1 . Flow streamlines and pressure contours for a simulation with IMF Bz>0. The 
pressure contours were placed in the noon-midnight and equatorial planes. The flow 
streamlines were started in the solar wind. The bottom panel shows an enlargement of the 
region close to the Earth. Additional streamlines have been added near the Sun-Earth 
line. 
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Flow Streamlines and Magnetic Field Lines 

(Bz>0, B y =0) 



Figure 2. Magnetic field lines calculated along two flow streamlines from the simulation 
in Figure 1. The colored numbers indicate the locations along the streamlines discussed in 
the text. 
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Field Lines Along a Dusk LLBL Ion Trajectory 

(B z >0; B y >0) 



Figure 3. Field lines (purple) along the trajectory of a solar wind ion (white). The 
trajectory was calculated in the electric and magnetic fields from a simulation with IMF 
Bz>0 and By>0. At this time B z =4.8nT and Bv = 1.3nT. Pressure contours have been 
plotted in the equatorial plane. They have been made semi-transparent. The orange line is 
a flow streamline started from the same position as the ion trajectory. 
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Field Lines Along a Dusk LLBL Ion Trajectory 

(B z >0; B y >0) 



Figure 4. Two views of the ion trajectory in Figure 3. The view in the top panel is from 
the dusk side while the view in the bottom panel is from the dawn side. Magnetic field 
lines have been calculated along the trajectory. The color coding on the field lines gives 
the magnitude of the magnetic field. 
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Figure 5. Projections of the ion trajectory in Figure 3 onto three planes (Z=0, Y=0, X=0). 
The trajectory is color coded with the particle energy. 
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Field Lines Along an Ion Trajectory 

(B >0; B >0) 



-10 



Figure 6 . Field lines (purple) along the trajectory of a solar wind ion (white). The 
trajectory was calculated in the electric and magnetic fields from the same simulation as 
Figure 3. Pressure contours have been plotted in the equatorial plane. They have been 
made semi-transparent. The orange line is a flow streamline started from the same 
position as the ion trajectory. 
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Ion Energy Along Trajectory 
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Figure 7. Projections of the ion trajectory in Figure 6 onto three planes (Z=0, Y=0, X=0) 
The trajectory is color coded with the particle energy. 






Ion Entry Locations 
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Figure 8. The distribution of ion entry locations on the dayside magnetopause. The top 
panel shows the entry locations projected onto the YZ plane while the bottom panel 
shows the locations projected onto the XY plane. The MHD simulation had IMF B z >0. 
(From Richard et at., 1994], 
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Flow Streamlines and Pressure Contours 



Figure 9. Flow streamlines and pressure contours for a simulation with IMF B z <0. The 
pressure contours were placed in the noon-midnight and equatorial planes. The flow 
streamlines were started in the solar wind and in the middle and inner magnetosphere. 
The bottom panel shows an enlargement of the region close to the Earth. Additional 
streamlines have been added near the Sun-Earth line. 
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Flow Streamlines and Magnetic Field Lines 

(B z <0, B y =0) 



Figure 10. Magnetic field lines calculated along three flow streamlines from the 
simulation in Figure 9. The colored numbers indicate the locations along the streamlines 
discussed in the text. 
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Magnetic Field Lines Along an Ion Trajectory 

(B z <0, By=0) 



Figure 11. Field lines (purple) along the trajectory of a solar wind ion (white). The 
trajectory was calculated in the electric and magnetic fields from a simulation with the 
IMF Bz<0. Pressure contours have been plotted in the equatorial plane. They have been 
made semi-transparent. The orange line is a flow streamline started from the same 
position as the ion trajectory. 
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